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We  calculate  the  energy  decay  rate  of  Josephson  qubits  and  superconducting  resonators  from 
nonequilibrium  quasiparticles.  The  decay  rates  from  experiments  are  shown  to  be  consistent  with 
predictions  based  on  a  prior  measurement  of  the  quasiparticle  density  nqp  =  10/yU.m3,  which  suggests 
that  nonequilibrium  quasiparticles  are  an  important  decoherence  source  for  Josephson  qubits.  Calculations 
of  the  energy-decay  and  diffusion  of  quasiparticles  also  indicate  that  prior  engineered  gap  and  trap 
structures,  which  reduce  the  density  of  quasiparticles,  should  be  redesigned  to  improve  their  efficacy.  This 
model  also  explains  a  striking  feature  in  Josephson  qubits  and  resonators — a  small  reduction  in  decay  rate 
with  increasing  temperature. 

DOI:  10. 1 103/PhysRevLett.  103.097002  PACS  numbers:  74.81. Fa,  03.65.Yz,  74.25.Nf,  74.50.+r 


Superconducting  integrated  circuits  are  a  leading  candi¬ 
date  for  scalable  quantum-information  processing  [1],  with 
recent  experiments  showing  accurate  state  control  of 
coupled  circuits  [2],  Josephson  qubits  are  relatively 
straightforward  to  couple  because  their  relatively  large 
size  enables  quantum  interactions  via  simple  wiring. 
However,  these  same  wires  also  make  the  qubit  susceptible 
to  various  decoherence  mechanisms,  and  thus  a  detailed 
understanding  of  all  loss  mechanisms  is  critically  needed 
[3].  Nonequilibrium  quasiparticles  are  known  to  be  an 
important  decoherence  mechanism  for  charge  qubits,  since 
quasiparticles  hopping  on  and  off  the  qubit  island  produce 
an  unpredictable  change  in  state  [4,5].  How  do  quasipar¬ 
ticles  affect  the  coherence  of  charge-insensitive  qubits, 
such  as  phase  [6],  flux,  and  transmon  [7]  devices? 

In  this  Letter,  we  present  a  simple  and  general  model 
that  describes  how  nonequilibrium  quasiparticles  produce 
energy  decay  in  superconducting  qubits  and  resonators.  We 
relate  predictions  to  the  quasiparticle  density  nqp:  based  on 
its  prior  measurement  [8],  we  calculate  an  energy  decay 
rate  T i  that  could  possibly  explain  current  limits.  Although 
other  decay  mechanisms,  such  as  surface  dielectric  loss 
[9],  are  important,  experimentalists  need  a  detailed  theory 
of  all  mechanisms  in  order  to  properly  design  experiments 
that  separate  out  and  measure  energy  decay,  and  to  ulti¬ 
mately  improve  qubit  performance.  We  conclude  that  qua¬ 
siparticles  must  be  considered  in  present  experiments. 

Our  model  is  described  in  four  parts.  First,  we  calculate 
the  energy  decay  rate  for  a  Josephson  qubit,  using  a  gen¬ 
eralization  of  the  environmental  P(E)  theory  [10].  Since 
the  prediction  is  based  on  tunneling  of  nonequilibrium 
quasiparticles,  we  next  calculate  their  density  versus  en¬ 
ergy  from  standard  electron-phonon  rates,  showing  the 
injection  energy  is  unimportant.  Third,  we  estimate  quasi¬ 
particle  diffusion  lengths  and  comment  on  the  design  of 
effective  quasiparticle  traps.  Finally,  we  consider  and  dis¬ 


cuss  possible  sources  of  nonequilibrium  quasiparticles, 
using  energy  estimates. 

No  equilibrium  quasiparticles  are  expected  in  a  super¬ 
conductor  well  below  the  critical  temperature  T  «  T(, 
because  of  the  exponential  suppression  of  excitations 
from  the  superconducting  gap  A.  However,  an  experiment 
on  Cooper-pair  boxes  has  measured  a  significant  quasipar¬ 
ticle  density  nqp  ~  10//rm3,  arising  from  an  unknown 
source  [8].  Once  excited,  quasiparticles  can  readily  tunnel 
through  a  Josephson  junction. 

Charge  associated  with  quasiparticles  allows  a  tunneling 
event  to  couple  energy  between  the  quasiparticle  and  the 
qubit,  as  described  by  the  environmental  P(E)  theory  [10]. 
When  a  qubit  state  of  energy  £j0  decays,  its  energy  is 
added  to  the  final  quasiparticle  state.  Energy  is  transferred, 
however,  in  only  a  fraction  of  the  tunneling  events  for  the 
case  of  large  junction  capacitance  C  that  we  consider  here. 
For  the  simple  case  of  the  transfer  of  charge  q,  the  proba¬ 
bility  for  qubit  decay  from  a  single  tunneling  event  is  p  — 
( q2/2C)/El0 .  We  extend  the  P(E)  theory  to  account  for 
quasiparticle  tunneling  having  both  electronlike  (q  =  —e) 
and  holelike  (q  =  e)  charge  transfer,  which  gives  a  qubit 
decay  probability  [11,12] 

e2/2  C,  „ 

P  ~  —p - (U/Ur  +  V,VrY,  (1) 

£io 

where  u  and  v  are  the  (energy  dependent)  BCS  amplitude 
factors  for  the  left  (/)  and  right  (r)  junction  electrodes. 

For  a  tunnel  junction  with  normal  state  resistance  RT ,  the 
total  qubit  decay  rate  is  obtained  by  summing  over  all 
possible  quasiparticles  and  accounting  for  filled  and  empty 
states.  For  the  tunneling  of  quasiparticles  from  left  to  right, 
we  integrate  over  all  possible  initial  energies  E  to  obtain  a 
rate 
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Fi=r^  JA  PpiPrf,[1~fr]dE  (2) 

1  f°°  dE  E(E  +  El0)  +  A2 

RtC  J a  Ew  Ve2  -  A 2V(£  +  E io)2  -  A2 
X  /;(E)[1  -  fr(E  +  El0)l  (3) 

where  piir(E)  =  E/{E2  —  A 2)1//2  is  the  normalized  quasi¬ 
particle  density  of  states,  and  ftr  is  the  nonequilibrium 
occupation  probability.  The  factors  ur,  vr,  pr,  and  fr  are 
computed  for  the  final  quasiparticle  state  with  energy  E  + 
El0,  and  the  coherence  factor  is  opposite  in  sign  as  for 
charging  effects  [8].  The  total  decay  rate  T ,  is  obtained  by 
also  summing  the  quasiparticle  tunneling  rate  from  right  to 
left  T  [ .  We  note  the  similarity  of  this  formula  to  the  Mattis- 
Bardeen  theory  for  metallic  conductors  [13,14], 

The  total  density  of  nonequilibrium  quasiparticles  is 
given  by  nqp  =  2 D(EF)  f  p(E)  f(E)dE,  where  D(EF)  is 
the  density  of  electron  states  at  the  Fermi  energy.  The 
integral  for  the  decay  rate  can  be  solved  for  the  physically 
relevant  case  of  contributions  from  f(E )  only  near  the  gap, 
as  motivated  below.  For  small  nqp,  the  total  decay  rate  [15] 
is 


72  /  A  \3/2  nqp 

RtC\eJ  D(Ef)  A' 


(4) 


The  decay  rate  scales  as  the  fractional  quasiparticle  occu¬ 
pation,  which  is  normalized  by  the  total  density  of  Cooper 
pairs  D(EF) A. 

We  can  immediately  check  the  relevance  of  this  predic¬ 
tion  by  plugging  in  typical  parameters  for  qubits.  In 
Ref.  [8],  a  quasiparticle  density  of  1 0/ pm2  was  experi¬ 
mentally  measured  in  charge  qubits.  For  the  aluminum 
superconductor  that  is  typically  used,  we  have  A/k  = 
2.1  KandD(£y)A  =  2.8  X  106//rm3.  Typical  phase  qubit 
parameters  are  El0/h  =  6  GHz,  R,  =  200  11,  and  C  = 
1  pF,  which  yields  from  Eq.  (4)  a  qubit  decay  rate  T l  = 
1/(2. 1  p s).  This  rate  is  reasonably  close  to  what  is  mea¬ 
sured  for  our  qubit  [16]  T {  ^  550  ns,  especially  consider¬ 
ing  Ref.  [8]  possibly  had  lower  nqp  because  of  quasiparticle 
traps. 

Nonequilibrium  quasiparticles  also  produce  energy  de¬ 
cay  in  resonators  from  loss  in  superconducting  wires. 
Following  the  calculations  in  Ref.  [17],  the  resonator  qual¬ 
ity  factor  Q  can  be  expressed  as  a  function  of  the  quasi¬ 
particle  density  [12].  Here,  the  bulk  conductivity  of  the 
superconductor  is  given  by  the  Mattis-Bardeen  theory, 
from  which  the  surface  impedance  can  be  expressed  in 
various  limits.  Geometric  factors  expressing  the  nonuni¬ 
form  current  density  in  the  coplanar  transmission  line  are 
also  included,  to  yield  in  the  thick-hlm  limit 

A(AY,2A!sl  (5) 

Q  s  gm  r  tt  V&v  D(Ep) A 


Here,  A  ^  (77/r(l<T„  A/ft)-  1/2  is  the  penetration  depth  cor¬ 
responding  to  the  surface  kinetic  inductance  /r0A,  crn  is  the 
normal  state  conductivity,  s  is  the  width  of  the  center 
conductor,  g  =  1.0  and  gm  =0.31  are  geometrical  factors 
expressing  the  effect  of  nonuniform  current  density  on  the 
kinetic  and  magnetic  inductance,  respectively,  and  y  = 
1/2  is  a  factor  concerning  the  interface.  Using  s  = 
3  p, m,  A  =  50  nm  for  Al,  a  resonator  frequency  of 
6  GHz,  and  quasiparticle  parameters  described  previously, 
we  compute  Q  —  107,  which  corresponds  to  a  decay  time  of 
300  ps.  This  time  is  significantly  longer  than  found  above 
for  qubits,  and  agrees  with  experimental  observation  that 
decay  times  are  longer  for  resonators  [18]. 

A  more  detailed  understanding  of  nonequilibrium  qua¬ 
siparticles  can  be  obtained  by  calculating  their  occupation 
probability  f(E)  using  electron-phonon  scattering  rates 
[19].  We  first  consider  a  bulk  superconductor.  The  lifetime 
of  a  quasiparticle  at  energy  e  to  scatter  to  any  energy  e'  by 
emitting  or  absorbing  a  phonon  of  energy  e  —  e'  is 


T'  ,  =  — 

^  To  4  A 


L 


n  2 


(>-£) 


X 


( kTc )3  Ve'2  -  A2 

i  -m 


exp[-(e  -  e')/kT  ]  - 


(6) 


where  r0  ~  400  ns  is  the  measured  characteristic  electron- 
phonon  time  for  Al.  The  lifetime  of  a  quasiparticle  state 
with  energy  e  to  recombine  with  another  quasiparticle  of 
any  energy  e'  is 


To  J  A  (kTcf  7  e'2  -  A2  \  ee'J 

1 

X  I  exp[  — (e  +  e')/kTp\  -  1|‘ 


Note  the  similarity  in  this  formula  to  Eq.  (6),  but  here  the 
energy  of  the  emitted  phonon  is  greater  than  2A.  With  the 
term  f(e '),  this  rate  is  proportional  to  the  density  of  qua¬ 
siparticles,  implying  that  the  recombination  rate  slows 
down  once  the  density  of  quasiparticles  becomes  small. 

Equations  (6)  and  (7)  are  used  to  solve  for  the  steady- 
state  occupation  probability  f(E)  of  quasiparticles  for  a 
given  (nonequilibrium)  injection  rate  rqp.  With  the  recom¬ 
bination  rate  being  proportional  to  the  square  of  the  density 
of  quasiparticles,  we  solve  this  nonlinear  differential  equa¬ 
tion  numerically.  Because  the  quasiparticle  relaxation  rate 
at  high  energy  is  dominated  by  electron-phonon  scattering, 
which  is  number  conserving  for  quasiparticles,  the  exact 
energies  at  which  the  quasiparticles  are  injected  are  unim¬ 
portant  for  the  computed  low-energy  distribution.  For  sim¬ 
plicity,  we  inject  quasiparticles  with  energies  between 
2.8A  and  3A. 

We  plot  the  quasiparticle  occupation  probability  versus 
energy  in  Fig.  1  for  an  injection  rate  rqp/2Z9(Ef)A  = 
1.7  X  l(L  l0/rn  that  gives  a  total  quasiparticle  density 
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FIG.  1.  Plot  of  quasiparticle  occupation  versus  energy  for  3 
phonon  temperatures  Tp  =  140,  70,  and  0  mK  (top  to  bottom). 
An  injection  rate  was  chosen  to  match  the  experimental  quasi¬ 
particle  density  10//zm3,  equivalent  to  n^/lDiEpiA  = 
SfpdE  =  1.8  X  10  6.  The  state  occupation  at  energies  shown 
does  not  depend  on  the  injection  energies.  The  gray  line  shows 
the  thermal  occupation  fT  at  140  mK. 

nj2D{EF)A  =  1.8  X  10  6,  approximately  the  same  as 
for  Ref.  [8],  At  phonon  temperatures  T  S  170  mK,  a  ther¬ 
mal  occupation  fT  of  quasiparticles  dominates.  At  T  ~ 
140  mK,  an  exponentially  decaying  occupation  is  calcu¬ 
lated,  but  at  an  elevated  occupation  due  to  the  nonequilib¬ 
rium  quasiparticles.  For  temperatures  below  about  70  mK, 
the  curves  become  independent  of  temperature.  This  cal¬ 
culation  is  consistent  with  experiments  on  Cooper-pair-box 
experiments  [8,20],  where  nonequilibrium  quasiparticle 
tunneling  rates  became  independent  of  refrigerator  tem¬ 
perature  below  about  70  mK. 

The  total  quasiparticle  density  at  T  =  0  is  numerically 
foundtobenqp/2D(£'/r)A  =  0. 14[r0rqp/2D(£'F)A]1/2.  We 
find  the  computed  decay  rate  is  well  approximated  by 
Eq.  (4). 

The  calculated  temperature  dependence  of  the  decay 
rate  is  plotted  in  Fig.  2(a),  which  shows  negligible  change 
at  low  temperatures  T  £  40  mK,  and  then  a  small  reduc¬ 
tion  in  rate  up  to  T„  =  120  mK.  This  unusual  prediction  of 
an  improvement  in  decay  rate  with  increasing  temperature 
results  from  nonequilibrium  quasiparticles  having,  on  av¬ 
erage,  larger  energy  and  thus  smaller  contributions  from 
density  of  states  and  coherence  factors.  At  higher  tempera¬ 
tures,  the  thermal  generation  of  quasiparticles  rapidly  in¬ 
creases  the  decay  rate.  This  temperature  dependence  well 
explains  data  taken  on  superconducting  resonators  [21]. 
Measurements  of  T^r)  for  a  phase  qubit  [16]  are  shown  in 
Fig.  2(b),  where  a  reduction  and  then  an  increase  in  decay 
rate  is  in  qualitative  agreement  with  predictions  [22],  The 
simplicity  of  our  model,  which  does  not  account  for  the 
small  occupation  of  the  qubit  ground  state,  diffusion  of 


0  0.05  0.1  0.15  0.2 

Temperature  Tp  [K] 

FIG.  2.  (a)  Plot  of  calculated  qubit  energy  relaxation  rate 

versus  phonon  temperature  Tp  for  a  quasiparticle  density  of 
nqp/2Z)(£'f)A  =  1.8  X  10-6  and  El0/h  =  6  GHz.  This  shape 
of  the  curve  is  relatively  insensitive  to  changes  in  Ei0  and  n^. 
Note  the  small  decrease  in  decay  rate  with  increasing  tempera¬ 
ture  up  to  120  mK,  and  then  a  rapid  increase  due  to  thermally 
generated  quasiparticles,  (b)  Experimental  measurement  of  re¬ 
laxation  rate  F|  for  a  phase  qubit.  The  relaxation  decreases,  and 
then  increases  with  temperature,  in  qualitative  agreement  with 
theory.  The  slow  change  at  T p  S  140  mK  is  probably  due  to 
nonidealities,  such  as  the  small  occupation  of  the  ground  state 
(greater  than  90%  for  Tp  <  100  mK,  but  only  77%  at  Tp  = 
130  mK  and  55%  at  T p  =  160  mK). 


quasiparticles,  or  gap  inhomogeneities,  may  explain  the 
differences. 

Previous  experiments  have  used  engineered  gap  [23]  and 
trap  [4,24]  structures  to  reduce  the  density  of  quasipar¬ 
ticles.  Are  they  large  enough  to  be  completely  effective? 
The  quasiparticle  diffusion  constant  for  aluminum  is  D  = 
60uqp  cm2/s,  where  uqp  =  (1  —  A 2 / £i2) 1  /2  is  the  normal¬ 
ized  velocity,  which  depends  on  energy  and  approaches 
zero  at  the  gap.  Quasiparticles  diffuse  a  distance  (Dr)1//2  in 
time  r;  we  plot  this  length  in  Fig.  3  for  both  the  electron- 
phonon  and  electron-electron  decay  times.  We  estimate 
that  the  characteristic  distance  for  quasiparticles  to  equili¬ 
brate  are  of  millimeter  lengths,  much  larger  than  the  size  of 
qubit  devices.  In  normal  metal  traps,  quasiparticles  lose 
energy  at  much  shorter  length  scales,  of  order  30  /zm. 
Finally,  quasiparticle  energy  is  ultimately  removed  from 
the  superconductor  by  recombination  and  emission  of  a 
phonon  with  energy  S  2A.  Because  phonons  are  long- 
lived  at  low  temperature,  especially  for  crystalline  sub¬ 
strates,  they  may  ballistically  travel  across  the  chip  and 
have  their  energy  redeposited  as  quasiparticles  anywhere  in 
the  superconductor.  All  of  these  estimates  imply  that  a 
conservative  design  should  have  quasiparticle  trapping 
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FIG.  3.  Plot  of  average  quasiparticle  diffusion  length  versus 
energy  for  both  the  electron-phonon  1  / 1 (e-ph)  and  electron 
recombination  l/L  ( e-e )  decay  times.  The  e-e  interaction 
dominates  for  E  S  1.1  A.  Parameters  are  D  =  60uqp  cm2/s 
(Al),  Tp  =  0,  and  =  10//rm3. 

structures  of  millimeter  or  greater  size,  not  just  local 
structures  placed  around  the  tunnel  junctions. 

Although  the  density  of  quasiparticles  depends  on  ex¬ 
perimental  details  of  chip  design,  it  is  possible  to  make 
rough  power  estimates  to  help  understand  the  generation 
mechanism.  The  model  presented  here  found  a  quasi¬ 
particle  injection  rate  rqp  ~  2.4  X  103/s  /xm  !.  Assuming 
a  total  superconductor  volume  of  0. 1  /im  X  10  mm2  and 
an  energy  A  per  quasiparticle,  the  total  power  to  the  chip  is 
0.06  pW. 

The  power  load  from  various  sources  can  be  compared 
with  this  number.  Cosmic  rays  have  a  flux  of  ~0.6/cm2  s 
and  deposit  an  energy  of  ~1  MeV/mm,  yielding  for  a 
50  mm2  chip  of  thickness  0.5  mm  a  power  0.02  pW. 
Background  radioactivity  is  typically  of  order  that  coming 
from  cosmic  rays.  If  coaxial  lines  allow  ~4  K  thermal 
radiation  to  be  transmitted  with  a  bandwidth  100  GHz,  a 
power  ~5  pW  is  found.  Blackbody  radiation  from  1  K 
gives  about  60  pW  of  power  load  to  the  outside  of  the  chip 
mount.  Materials  also  slowly  release  energy  at  low  tem¬ 
peratures,  with  amorphous  SiOi,  poly  crystalline  Al,  and 
Teflon  showing  heat  release  of  —100  pW/g,  —20  pW/g, 
and  —2  pW/g,  respectively,  that  slowly  decays  with  time 
[25].  From  these  estimates,  nonequilibrium  quasiparticles 
are  clearly  plausible. 

In  conclusion,  we  have  developed  a  model  that  predicts 
qubit  energy  decay  from  the  density  of  nonequilibrium 
quasiparticles.  Predictions  of  the  model  are  in  reasonable 
agreement  based  on  a  prior  experimental  measurement 
and,  at  a  minimum,  show  that  this  decoherence  mechanism 
should  be  carefully  considered.  The  quasiparticle  energy 
distribution  calculated  from  simple  electron-phonon  inter¬ 
action  qualitatively  agrees  with  measurements,  and  we 
argue  that  engineered  gap  and  trap  structures  should  be 
of  large  (millimeter)  size.  With  a  good  understanding  of  the 
basic  physics,  we  believe  experiments  can  be  performed  to 
more  carefully  test  for  nonequilibrium  quasiparticles,  and 


to  ultimately  improve  the  coherence  of  superconducting 
qubits. 
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